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ABSTRACT 

The instability of hypersonic boundary- layer flows over flat plates is considered. The 
viscosity of the fluid is taken to be governed by Sutherland’s law, which gives a much more 
accurate representation of the temperature dependence of fluid viscosity at hypersonic speeds 
than Chapman’s approximate linear law; although at lower speeds the temperature variation 
of the mean state is less pronounced so that the Chapman law can be used with some con- 
fidence. Attention is focussed on the so-called “vorticity mode of instability of the viscous 
hypersonic boundary layer. This is thought to be the fastest growing inmscid disturbance 
at hypersonic speeds; it is also believed to have an asymptotically larger growth rate than 
any viscous or centrifugal instability. As a starting point we investigate the instability of 
the hypersonic boundary layer which exists far downstream from the leading edge of e 
plate/ In this regime the shock that is attached to the leading edge of the plate plays no 
role, so that the basic boundary layer is non-interactive. It is shown that the vorticity mode 
of instability^ this flow operates on a significantly different lengthscale than that obtained 
if a Chapmanviscosity law is assumed (see Smith and Brown, 1989). In particular, we find 
that the growth rate predicted by a linear viscosity law overestimates the size of the growth 
rate bv 0(M 2 ). Next, the development of the vorticity mode as the wavenumber decreases is 
described, and it is shown that acoustic modes emerge when the wavenumber has decreased 
from it’s 0(1) initial value to 0(M"i). Finally, the inviscid instability of the boundary 
layer near the leading edge in the interaction zone is discussed and particular attention is 
focussed on the strong interaction region which occurs sufficiently close to the leading edge. 
We find that the vorticity mode in this regime is again unstable, and that it is concentrated 
in the transition layer at the edge of the boundary layer where the temperature adjusts 
from its large, 0(M 2 ), value in the viscous boundary layer, to its 0(1) free stream value. 
The existence of the shock indirectly, but significantly, influences the instability problem y 
modifying the basic flow structure in this layer. ■ • 


1 Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 

Our concern is with the inviscid mode of instability of hypersonic boundary-layer flows 
over flat plates. In the first instance we will consider a regime where there are no shocks 
present, and then we will show how the instability problem is significantly modified by 
the presence of a shock in the flow field. The motivation for this and related work on 
hypersonic boundary-layer instability theory is the renewed interest in hypersonic flight 
which has been stimulated by plans to build a successor to the Space Shuttle. A primary 
concern with such a vehicle is the question of where transition will occur over a wide 
range of Mach numbers and whether it can be controlled. At the largest relevant Mach 
numbers, say Mach 20-25, the extremely high temperatures associated with the flow would 
destroy the vehicle unless it were cooled, so that it is of interest to know the effect of 
the wall temperature on the instability properties of the flow. The purpose of this paper 
is to determine the inviscid instability characteristics of physically realistic hypersonic 
boundary-layer flows. We note here in passing that there is a simple generalization of 
Rayleigh’s (incompressible) inflection point theorem to compressible flows (Lees & Lin, 
1946), and that many compressible boundary layers turn out to be inviscidly unstable even 
though their incompressible counterparts are stable. This is a significant result because 
the growth rates of inviscid disturbances are generally much larger than those of viscous 
or centrifugal instabilities so that they are the likely cause of transition to turbulence 
in most situations. The modes which we discuss in this paper can be referred to as 
generalized inflection point modes because, when neutral, their phase speed is equal to the 
fluid velocity at the generalized inflection point. Furthermore the eigenfunctions of the 
modes are localized around that point. 

For convenience we will concentrate on high-Reynolds-number flow past a flat plate, 
although many aspects of our analysis are applicable to other boundary-layer flows (e.g. 
flow past a wedge). The underlying steady flows that we study depend upon the Mach 
number, the Prandtl number and the choice of the viscosity law; the complications arising 
from real gas effects are not investigated. Throughout we assume that the fluid viscosity 
is Newtonian and is adequately described by Sutherland’s formula. We wiU also take 
the Prandtl number to be one, noting that it is relatively straightforward to relax this 
restriction. Moreover the relaxation of the latter restriction does not significantly alter the 
qualitative features of the results we present here. 
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Reshotko (1976) and Mack (1987) have reviewed earlier work on the linear instability 
of high-Reynolds-number compressible flows. Many of these studies are based on the Orr- 
Sommerfeld equation; for a critique of the mathematical rigor of this approach see Smith 
(1979, 1989). Here we examine the linear stability of high-Reynolds-number flows by means 
of formal asymptotic expansions; for example Smith (1989) has applied triple-deck theory 
to the lower-branch viscous Tollmien-Schlichting modes of compressible boundary layers. 
Seddougui, Bowles & Smith (1989) have extended this theory to include the effects of 
severe wall cooling, while Cowley & Hall (1988) have shown how such modes can interact 
with a shock at large Mach number. However, the viscous modes have relatively small 
growth rates, and our main concern here will be with inviscid modes. The nature of the 
asymptotic expansion procedure in these investigations clearly depends on the nature of 
the mode of instability. In fact the third type of instability responsible for boundary- 
layer transition, the Gortler vortex mode, develops an asymptotic structure at high Mach 
numbers closely related to that found here; see Hall and Fu (1989). 

When a quasi-parallel approximation is formally justifiable because the Reynolds num- 
ber is large, inviscid modes satisfy the compressible Rayleigh equation. Numerical solutions 
to this equation have been reported by, inter alia, Mack (1984, 1987) for boundary-layer 
flows, Grosch & Jackson (1989) for shear flows, and Papageorgiou (1989) for wake flows. 
For fluids satisfying a Chapman viscosity law, high-Mach-number asymptotic solutions to 
this equation for the so-called “acoustic” boundary-layer modes have been obtained by 
Cowley &: Hall (1989), while Smith & Brown (1989) have identified the asymptotic form 
of the “vorticity” mode. Goldstein & Balsa (1989) have given an asymptotic solution for 
high-Mach-number, shear-layer modes of instability for a Chapman fluid. Though the basic 
states investigated by the two latter pairs of authors are different, they found essentailly 
the same most unstable eigenvalue because it corresponds to a disturbance trapped in a 
thin layer where the overall features of the basic state are unimportant. However Goldstein 
and Balsa did not spot the exact solution of the vorticity mode equation found by Smii;h 
and Brown. 

Both in the above mentioned boundary-layer analyses, and the hypersonic Gortler 
vortex instability analysis of Hall & Fu (1989), one of the key asymptotic regions for the 
case of a Chapman viscosity law is a logarithmically thin layer which develops due to the 
exponential decay of the underlying steady temperature field away from the wall. However, 
Chapman’s viscosity law is not exact, and was introduced as a useful interpolation law 
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which greatly simplified steady boundary-layer calculations (e.g. see Stewartson 1964). At 
the large temperatures typical in hypersonic flows, it differs significantly from the more 
precise Sutherland’s formula. In fact Chapman’s law is simply a linear approximation to 
the viscosity-temperature dependence of the fluid; it is therefore of questionable validity 
in the hypersonic limit. At high Mach numbers the steady temperature field in a fluid 
satisfying Sutherland’s formula initially decays algebraically away from the wall, before 
reverting to exponential decay in an asymptotic region “far” from the wall (e.g. Freeman 
& Lam 1959). The effect of this algebraic decay is to change significantly the scalings in the 
transition region; in particular the asymptotic expansions proceed in inverse powers of M 
rather than \/log(M). Moreover the wavelength of the most unstable disturbance varies 
by a factor of \/log(M) in the two cases. We note that a similar difference in scalings is 
evident in the interaction region of steady hypersonic flow past a flat plate. In that case 
Lee & Cheng (1969) have shown that the shock-heating transition layer is logarithmically 
thin for Chapman’s viscosity law, whereas for a power-law viscosity formula, and hence 
for Sutherland’s formula, the scaling for the transition layer is algebraic (Bush 1966). 

The flows which we consider here are appropriate to different stages of hypersonic flow 
past a semi-infinite flat plate. In the first instance we shall consider the instability of a 
non-interactive flow which is appropriate to large distances downstream of the leading edge 
of the plate. Here the attached shock at the leading edge has no effect on the flow field and 
the basic state is the Sutherland law counterpart of that discussed by Smith and Brown 
(1989). This basic state, and the Rayleigh equation which governs the inviscid instability, 
are discussed in Section 2. The dispersion relationship associated with this equation is 
then solved in Section 3. We shall discuss the growth rate of the mode over the whole 
range of unstable wavenumbers and discuss how the disturbance is related to the acoustic 
mode at small wavenumbers. 

Then in Section 4 we go on to discuss the basic state in the so-called “strong interaction 
region” further upstream. The description of the basic flow in this regime for a power law 
fluid is due to Luniev (1959) and Bush (1966). We shall discuss how the interactive system 
for the flow can be formulated when the Mach number is of order R*. However this 
system can only be solved numerically and we are not aware of any published results on 
that problem. Nevertheless we can of course still consider the instability of that state, 
and we derive the appropriate (quasi-parallel) stability equations in this regime. The 
strong hypersonic interaction limit then corresponds to taking the further limit of the 
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streamwise variable tending to zero. In that case a similarity solution for the basic fi.ow 
can be found, Bush (1966), and a re-scaled Rayleigh equation for the disturbance is found. 
The solution of that equation is discussed in Section 5. We could have instead considered 
the weak hypersonic limit further downstream where Bush and Cross (1967) have given an 
appropriate asymptotic description. We choose to concentrate on the strong interaction 
regime because it is, to a certain extent, simpler. Further, if the flow is unstable in this 
regime it is likely that growing disturbances will originate here. Finally in Section 6 we 
shall draw some conclusions. 
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2. NON-INTERACTIVE STEADY FLOWS 


We begin by considering the stability of steady hypersonic flow far downstream from 
any leading-edge interaction region; in particular, if L is the distance from the leading 
edge, and , a x , p x and £«, , are the velocity, sound speed, density and shear viscosity 
of the free stream flow, then we assume that the Reynolds number, 

R = Poo Uoo± (2.1a) 

Poo 

is larger than whatever power of the Mach number, 

M = 2 ^, (2-W) 


is necessary for interactive and/or non-parallel effects to be negligible (see below for a more 
precise restriction). We adopt a non-dimensionalisation based on coordinates Lx_ (where x 
is in the direction of flow and y is normal to plate), velocities UooH, time LtjlJ x , pressure 
p^Ulp, density temperature T^T, and shear and bulk viscosities PooP and p x p 

respectively, where the subscript oo denotes the value of the quantity in the free-stream. 
On the assumption that the fluid is a perfect gas with a constant ratio of specific heats 7 , 
the governing equations of the flow are 


where 


^ + v ■ (pu) = 0 , 

dt 

(2.2a) 

p^= = -Vj)+ -{pe) + V((p' ~^p)V -u)l 

Dt n * 

(2.26) 

= b-l)M^ +p >(,VT )+ ^ J #. 

(2.2c) 

7 M 2 p = pT, 

(2.2 d) 

__ 1 ( dui dUj \ 

2 ' dxj dXi ' ’ 

(2.3a) 

2 

$ = 2pg : g + (p 1 - -p)(V ■ u) 2 , 

(2.36) 


and Pr is the constant Prandtl number. For a shear viscosity obeying Sutherland’s law 


. _ 110.4 

where C w — - — for air temperatures measured in degrees Kelvin. In the numerical cal- 

1 oo 

culations discussed in the next section we took f x = 216.9. The boundary-layer equation 
can be recovered by first substituting 


-w 

J A 


pdy, v = R > V, 


(2.4a, 6) 


where the Dorodnitsyn-Howarth variable, C, is introduced for convenience, and then taking 
the limit R — » oo. 

For steady two-dimensional flow over a flat plate, a similarity solution to these equa- 
tions exists. With 


V ~ V / (1 + C)x’ PV ~ + 

rl>=y/(l + C)xf(ri), T = T(t]), p = p{rj) y 


(2.5 a) 


(2.56) 


the governing similarity equations are found to be 


pT = 1, 


(2.6a) 


5"" + (rT c f ") = 0 ' 


(2.66) 


l f r + ±(jj-r) + (^i)W 2 = 0 

2 7 Pr\T + Cj T+C 1 ’ 


(2.6c) 


subject to the boundary conditions 


/(0) = /'(O) = 0, 7'(oo) = r(oo) = 1, 


(2.6d) 


and T(0) = T w (fixed wall — temperature), or T'(0) = 0 (insulated wall). (2.6e) 

For simplicity we will focus attention on Pr = 1, and denote by T r the wall temperature 
when the boundary is insulated. Then, as is well known (e.g. Stewartson 1964), the energy 
equation can be integrated to yield 


T = 1 + ((r, - 1) + -(7 - 1 )M 2 (T„ + /'))(! - /'), 


where T w = T h T r and T r = 1 + ^(7 - 1 )M 2 . 


6 


nl WiH II 


The solution to (2.6) in the limit of large Mach number has been examined by Freeman 
& Lam (1959). They showed that two asymptotic regions develop, defined as where the 
coordinates r) and £ = M*rj, respectively, are order one. 


* = 0(D 


In this region we write / = M'*/o(*) + • • then using (2.7) it follows from (2.6 6) 


that 


with 


, r ,m-f A V = o, 

fofo + y 1 ) \(Ti + f>)Ui- 


(Tt + nHi-n- 

/o(o) = /6(°) = o, 


(2.8a) 


( 2 . 86 ) 


and 

, f. B . 24 -, + ... as? ^<x>. (2.8c) 

fo ~ { + f> + pT + 1)(7 _ ijjj ( Tl + l)( 7 - i)(4 

We note that as { -> oo, then fo -> 1 algebraicaUy. This is of course different than 
the corresponding result for a Chapman fluid, in that case the correction term to the 
free-stream speed is exponentially small. This difference is significant because it leads to 
inviscid instabilities which are unstable at very different wavelengths. 


77 = 0(1) 


Here, we write 


, , P , 91 

f = v + l^ + iP + 


then gi satisfies the equation 


m'{ + 2 


( i/i-ih-OPi + iKi „V 

\ 1 + C - $(1 - 1)(T. + l)pi 91 ) 


= 0, 


(2.9) 


(2.10a) 


subject to 

0 , — — w as 77 — > 0, and gi -* 0 as rj — ► oo. (2.106) 

91 (T, + 1)(7 - l)*? 3 

In Figure (2.1a) we show the function /q for the adiabatic case = 1 with 7 — 1.4, we 
note the algebraic approach of the function to the free-stream speed. In Figure (2.1b) 
we show the transition layer function g[ corresponding to the same case. We note the 
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exponential decay of this function for large tj. It is also worth pointing out that for a 
Chapman fluid the transition-layer equation corresponding to (2.10a) is linear, and it’s 
solution can be expressed in terms of the exponential function. It was this simplification 
that enabled Smith and Brown (1989) to spot the exact solution of the neutral vorticity 
mode in their study of the instability problem in this layer. We shall see below that for 
Sutherland’s law we are not able to find a similar exact solution of the stability equation. 

Sufficiently far downstream the quasi-parallel assumption is valid for inviscid instabil- 
ity modes. It is then appropriate to seek perturbations of the form 

u = f'(v) + • .. + Au(r?)exp (*y ( ] ~j~C) ( QX + & z ~ c *)) +•••» ( 2 . 11 a) 


P= ^M2 + -' + A ^ )eXP 


(ora: + fiz — ct ) 1 + . . . 


with similar expressions for the other flow quantites. If the disturbance amplitude, A, is 
sufficiently small the pressure perturbation p satisfies the compressible Rayleigh equation 


d 2 p _ 2/" dp 
drj 2 f — c dr) 


( a 2 + / 3 2 )T(T 


a 2 M 2 {f - c ) 2 
(a 2 + P 2 ) 


)P = 0 , 


(2.12a) 


where for convenience we have rescaled the quantities a, ft, c by dividing them by . The 
conditions that there is no normal velocity at the wall, and that the disturbance is confined 
to the boundary layer, can be expressed as 


p' = 0 on r) 


0 as 77 


(2.126) 


Equation (2.12a) and boundary conditions (2.126) specify a temporal-stability eigenrelation 
c = c(a,/5); alternatively from a spatial stability standpoint, the eigenrelation can be 
regarded as a = a(c,/ 3 ). 
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3. THE FAR DOWNSTREAM BEHAVIOUR OF THE INVISCID MODES 

In this section we discuss the asymptotic form of unstable solutions to (2.12a) for the 
region far downstream of the leading edge of the plate. In a previous investigation Cowley 
and Hall (1988), hereafter referred to as CH, studied the so-called acoustic modes of (2.12a) 
in this region on the assumption that the viscosity satisfies Chapman’s law. Simultaneously 
Smith and Brown (1989), hereafter referred to as SB, investigated the vorticity mode for 
Chapman’s law. The main difference between these two types of modes is that the acoustic 
modes have a ~ M~ 2 , whilst the vorticity mode, at least while it is close to neutral, has 
a ~ \/ 2 log M 2 . Moreover the vorticity mode is centred at the adjustment layer at the 
edge of the boundary layer whilst the acoustic one is concentrated in the main part of the 
boundary layer. 

However, as indicated above, at high Mach numbers the temperature variations in the 
boundary layer are large, so that a linear temperature- viscosity law is a bad approximation; 
Sutherland’s law should be used to give a better representation of the viscosity. It is then 
important to see how the asymptotic structures developed by CH and SB change. We shall 
see that there are significant differences. 

In the first instance we derive the asymptotic structure of the solution for a two- 
dimensional vorticity mode of (2.12a). We determine the neutral values of a and c for 
this mode, and find the limiting form of the mode when the further limit a — * 0 is taken. 
This limiting solution points to a sequence of distinguished asymptotic limits. Within the 
sequence of asymptotic limits mentioned above the scaling, a ~ M -3 ^ 2 , appropriate to 
an acoustic mode emerges; we therefore discuss the latter mode as a limiting case of the 
vorticity mode. 

a = 0(1) 

Consider then the solution of (2.12a) which has the eigenfunction trapped in the 
temperature adjustment layer at the edge of the boundary layer. We seek a solution which 
has /? = 0 and a = 0(1), so that the wavelength of the vorticity mode is comparable with 
the depth of the adjustment layer in terms of the Dorodnitsyn-Howarth variable. From 
(2.7) and (2.9), the velocity field, tz, and temperature field, T, of the underlying steady 
flow expand as 

20(77) 

U ~ 1+ (7-l)(T 6 + l)M 2 + '"’ (3.1a, b) 

T = 1 — G + . . . , 
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where 77 = y/2 Tj and the function G = ^(Ji, + 1)(7 — 1)^^ satisfies a rescaled version of 
(2.10a), and has the asymptotic behaviour 


G -~^ + 77 + (2C + 1) + . . . as T) 

G — ► 0 exponentially as Tj — * 00. 


0, 


Here B is a constant to be calculated numerically. 
Next we expand a, c and p in the forms 


a = — pd + . . . , 
v/2 

c = 1 + 


( 7 - i)(r, + 1 )m 2 

p = p+ ..., 


£ + • • - > 


(3.2a) 

(3.26) 


(3.3a, 6, c) 


where we have assumed that the disturbance moves downstream with the fluid speed in 
the adjustment layer. On substituting for u and T from (3.1), and using (3.3), we find 
that the zeroth order approximation to (2.12a) in the adjustment layer is, 


d 2 p 

dfj 2 


2 G' dp ^ 2 


- d z ( 1 - G) 2 p = 0 


G-c , 


(3.4) 


We note that three-dimensional disturbances satisfy (3.4) but with a replaced by 

n *7 t f f 

a + p so that in the present regime it is sufficient for us to consider only two-dimensional 
modes. Equation (3.4) is to be solved subject to p vanishing in the limits f) — * 0 and 
fj — * oo, so that the disturbance is confined to the adjustment layer. For fj >> 1 it follows 
from (3.4) that p decays like exp(— dfj), whilst for f) << 1 a WKB solution of (3.4) can be 
expressed in the form 

p ~ exp[— J Q(fj)]dfj, (3 5) 

where 

9d 2 Ba 

0 ~ t + — h — - — t = + 2 Ca + . . . as tj — » 0. (3 6) 

First, we restrict our attention to the neutral case, c is then real and can be evaluated 
by finding the fluid speed correct to order M ~ 2 at the generalized inflection point which 
is located where 

hi i- 2 S = o, 

U, T 
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i.e. where 


2 G' 



1 - G 


= 0 . 


A numerical solution to (2.10a) using a Runge-Kutta method shows that this occurs when 
7] ~ 1.971510, in which case 

c= -.993934 (3.7) 


The corresponding real value of a is obtained by integrating (3.4) from 7 ) = 0 to Tj = oo with 
an appropriate treatment at the generalised inflection point. Such a calculation predicts 
that the neutral value of a is 

a « 0.645065. (3.8) 

However, of greater significance are the unstable eigenmodes. Figure (3.1) illustrates 
the dependence of the growth rate, dc<, on the real wavenumber a . We see that the growth 
rate attains it’s maximum value of ~ 0.256853 at a ~ 0.143619; further it turns out that 
the acoustic modes have smaller growth rates (see below) so that this is the most unstable 
inviscid mode for a hypersonic boundary layer. In Figure (3.2) we show the eigenfunction 
of the vorticity mode equation at different values of the wavenumber. This figure indicates 
that as the wavenumber decreases the eigenfunction starts to expand out of the transition 
layer. We note here that c* as defined above is independent of T b and that the growth rate 
is obtained from (3.3b) by dividing ac { by ^(T* + 1 )(t - !)• It follows that wall cooling 
has a destabilizing effect on the vorticity mode. 

The structure of the growth-rate curves at small wavenumbers is of interest because 
for sufficiently small values of the wavenumber we expect that the vorticity mode will 
develop a structure similar to that of the acoustic mode. In the following discussion we 
will isolate the different significant regimes which occur in the small wavenumber limit. 
For simplicity we shall now concentrate on the adiabatic problem and take T b = 1. 

The key to understanding the subsequent regimes when a is related to inverse powers 
of the Mach number is to write down the small a asymptotic structure of (3.4). Figure 
(3.3) is a schematic illustration of the different regions in 7 ) space which emerge in the limit 
q. — > o. Also shown in this figure is the wall layer in which the temperature becomes of 
order M For the moment d is not considered to be sufficiently small for d to be 0(M ^) 
for some positive (j>\ it then turns out that the regions I and IV shown in figure (3.3) are 
quite passive. However, at sufficiently small values of d the wall layer structure of the 
basic state will enter the problem - see below. 
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After some careful numerical calculations at small values of d we deduced that c 
expands in the form 

ci C2 

+' — 1 - 

! A Q / *7 ' * • • • 


C “ d 4 / 7 + d 3 / 7 + ^ 3 ' 9 ^ 
Also, from (3. 2 b, 4) it follows that for f\ ~ d _1 the pressure decays like exp(— dfj). This 
suggests that the pressure in III should be expanded in the form 


P = 1 + d 7 pi + a~p2 + . . . + a'p6 + d(p7 — f?) + 

d*(p8 - Plf}) + 0t’(p9 ~ P2TJ) + A^(pio - P3#l)+ 
+ a~ Pii + a~ P\2 + d~Pi3 + d 2 Pi4 -f . . . , 


(3.10) 


where we have anticipated the form of several terms in this expansion in that pi,p2, . . . , Pio 
are taken as constants, whilst Pn,Pi 2 , etc. are functions of fj. At order d^ in this 
expansion we find that P\\ satisfies 


2 G' 
Cl 


Hi " ==- = o, 


so that after use of the exponential mathching condition it follows that 

2 f 00 

Pll = -— / GdTj - p4fj + pn, 
ci %- 


(3.11) 


where p\\ is another constant. P 12 and P 13 satisfy similar equations with forcing functions 
coming from the higher order terms in the expansion of the wavespeed. 

The order d 2 term in (3.10) is then found to satisfy 


Ha = (1 - GY 


2 G' , C2 „ ,£3 C2., C 4 2C2C3 cjL 

\P3 - —P2 —)P1 - — + — ^ , 

4 4 


c 1 ci ci cj Ci 

which may be integrated twice to give P 14 . However it is enough for our purposes to note 


that when fj —* 0, P 14 ~ fj 0 so that as fj — ► 0 the order d^ - and a 2 terms in (3.10) become 


comparable when 77 = 0(d?); thus, as anticipated earlier, the sublayer III is of depth d* 
Within this layer we define 


a 


and expand the pressure as 


p = 1 + a 'pi + . . . + dp 7 + a'P%(Z) + . . . , 


12 


where we have, by matching with the solution in II, again anticipated several terms in this 
expansion. At order a* we find that Ps satisfies 


Ps + 


72 

Z{c\Z 4 + 9) 
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We write the solution for Ps in the form 


^=(c 1+ |j) 2 (/ 


81dZ 


(diZ 4 + 9) 2 c\ 


-i). 


where the constant of integration has been chosen to satisfy Pg — ► — 1 as Z — ♦ oo in order 
to match with the solution in II. Finally we consider Region IV and define Tj = z. The 
zeroth order approximation to (3.4) here is 


8 81 

Pzz + ~Pz g P - 0, 

z z° 


(3.12) 


which must be solved subject to p — * 0 as z — ♦ 0. The appropriate solution then has p —* 
constant as z — > oo so a match with III can only be achieved if 



81 dZ 1 _ n 

(ciZ 4 + 9) a 


(3.13) 


It is assumed that c\ in (3.13) is complex so that after a little manipulation we obtain 


/ 8\/2 . ± 3i7r 

cl = - ( ^’ * M -r ] '• 


which corresponds to an unstable mode of (2.12a). 

a = Q(Af-») 

Next we consider the situation when a is so small that IV in Figure (3.3) merges with 
the wall layer of the basic state. At this stage the wall layer is of thickness M 3 so that 
a ~ M~ a . The wavenumber a is then written in the form 


where for the moment we shall not be precise about the size of the first correction term in 
this expansion. The zeroth order approximation to (2.12a) in the wall layer then becomes 


~n 

P 


2 of 


0 -1 


Go 


P ~ 


«o (7 - 1 ) : 


(1 + O 0 ) 2 (l - O 0 ) 2 (l - 


2(1 - 00) 


( 7 -l)(l+ Oo) 


— )p=0, (3.14) 
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where a dash denotes a derivative with respect to the wall layer variable £ = M» 17, and 
Co is the first term in the expansion of C in that layer. The above equation is again to be 
solved subject to p'(O) = 0. For large £ it has the asymptotic solutions 

n 

p ~ Nq = constant and p( ~ N\ = constant. 

For most values of a 0 the constant IVo is nonzero, and the structure in layers I, II, III 
survives intact. Thus for these values of ao the wavespeed c expands as 

d\ 

c = 1 - 77T + ••■, (3.15) 

and d\ will be complex so that the wave growth rate is of order M~ tt. 

However, equation (3.14) has a countable infinite set of eigenvalues for which the 
constant IVo = 0- In this case the eigenfunction takes on its greatest value in the wall 
layer, i.e. where the steady velocity field is adjusting to it’s free-stream value. These 
eigenvalues correspond to the acoustic inviscid modes and are the counterpart of those 
discussed by CH. A numerical solution of (3.14) yielded the sequence ao = 2.47, 7.17, 
12.19, 17.33, 22.54, 27.79, .... The first three eigenfunctions associated with this sequence 
of eigenfunctions are shown in Figure (3.4). Thus at a countable discrete set of points the 
acoustic modes emerge as almost neutral “vorticity” modes. In fact we could seek acoustic 
modes with a ~ 0(1), i.e. for a wavenumber much larger than currently assumed, for 
which the eigenfunction is concentrated in the £ = 0(1) wall layer. These modes have a 
fast variation in this layer which may be described by the WKB method. At certain values 
of M these eigenvalues coalesce with the neutral vorticity mode discussed earlier, and an 
analysis outlined in CH (see also SB) can be performed to describe the ‘splitting’ of the 
eigenvalues in this region, we do not pursue that calculation here but simply determine 
the structure of the vorticity mode in the vicinity of the acoustic mode. Thus the small a 
structure outlined above remains valid for a = 0(M~ % ) if the acoustic modes are avoided; 
it follows that c* continues to grow in the further limit aM » — ► 0, where no neutral acoustic 
modes exist. 

In order to examine the interaction of the vorticity mode with the acoustic modes in 
a little more detail, we expand the wavenumber in the form 

rv - Q 0 , <»1 , 

a - M 3/2 + M 25/U 
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whilst c still expands as in (3.15). Here we are denoting by «5 any eigenvalue of (3.14) 
corresponding to No = 0. In the wall layer p expands as 

p= M 2/7 Qo{1 + $[M)) + Qi + ( 3 - 16 ) 

where we have anticipated the possible occurrence of further eigenfunctions between Q o, Q l 
by inserting the factor (1 + $[M|). For our purposes it is not necessary for us to calculate 
$ here. If we substitute the above expansion into Rayleigh’s equation we see that Q o 
satisfies (3.14) with a = oq- The derivative of Qo vanishes at £ = 0 whilst for £ >> 1 

— 7 

Qo ~ qot 


where qo is a constant. At higher order in Rayleigh’s equation we find that Q\ satisfies 

2a' n . - 1 ) 2 . _ 2 W1 2(1 - Go) ^ 

Q i - — (! - u o)(l - L_ivi4^ ,yi 

Qq — 1 


a \ a 0 


(i - x)(i + «o) 
2(1 - G 0 )_ 


( 7 - 1) (1 - «o)(l - ( 7 _l)(l + ao ) 


)Q o- 


The solution of this equation which has Qi(0) = 0 is such that for large £ 


Ql ~ aoqoa*!- 


The constant g 0 depends on a' 0 and 7 , and can only be determined numerically. However it 
can be shown that ao is alternatively positive and negative at successive values of aft- The 
region III of Figure (3.3) becomes of depth 0(M~ 3/14 ) when a is 0(M~ 3 ' 2 ). The solution 
in this layer is calculated using the procedure outlined above; the only significant change 
is that in matching regions III and IV we must now account for the fact that Qo ~ T 7 
for large £. After some manipulation we find that the eigenrelation obtained when this 

matching has been carried out is 


j_ = i_ + l 

Ct\ d\ 


(3.17) 


where 

a\ = k = Oil, D = ( 9 ^> 4 ' 7 ’ A =(7rT)- (3.180,*, c. d) 

In Figure (3.5) we have shown solutions of (3.17) for -5 < a\ < 5. At both ends of this 
Figure Ji approaches e~ ^ which gives the required matching with the unstable vorticity 
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modes away from the acoustic modes. We further note that for small a| the growth rate 
tends to zero. More precisely from (3.17) we find that for small a\ 


di ~ + 


d\ 1/8 e - ir/ 4 


+ 


so that, dependent on on the sign of oj, the growth rate goes to zero like M~ or 

M~?(a i) 11/8 . 

Strictly we should now rescale the above expansions when aj becomes sufficiently 
small. We do not pursue that rescaled problem here since our main concern is with com- 
pleting a discussion of the structure of the vorticity mode at all lengthscales, in particular 
we wish to see if the vorticity mode connects with any other neutral states. We have 
already seen that Cj continues increasing as the wavenumber decreases through 0(il/ - ^^) 
values, apart from small neighborhoods of the acoustic modes where it decrease towards 
zero. Thus c, tends to infinity on the present scaling when ao tends to zero. Thus the 
mam significance of the a ~ M~ 3/2 regime is that it is at this stage where the acoustic 
modes emerge; however the small a structure, developed initially for the a = 0(1) scaling, 
survives this regime intact in the limit of o>q tending to zero. 


a = 

The next significant stage in the development of the vorticity mode occurs when the 
temperature in the upper layer, i.e. where t) ~ a -1 , becomes such that T ~ (o — c) 2 M 2 . 
The pressure eigenfunction in the upper layer then has its decay to zero modified. This 
situation occurs when 


c = 1 — 



(3.18) 


with 


q= mW + -"- (3.19) 

In this case the pressure in the upper layer decays like exp(— Z(1 — Jf)> This 
means that in the rescaled region corresponding to III of Figure (3.2), where we write 
Z = the pressure perturbation expands as 


and P 2 satisfies 


the equation 


p" 24 A 

2 Z[diZ 4 -3A} P2 


Z° 
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The solution of this equation which enables us to match with the behaviour discussed 
above is 

In the wall layer there are now no acoustic mode eigenfunctions which lead to decay at 
infinity so that matching leads to the eigenrelation 


2;2 


9 a 2 (7 - 1 yi 

(. hz 4 - 3 A) 2 


dZ — l\ 


1-d? 

d? 


)• 



9A 2 (7 - l) 2 
(JiZ 4 -3A) 2 


dZ + 


dl 


A more convenient form for the eigenrelation is found by writing 


8lV2 

1 ~ 3(3A) 1/4 (7 - 1) 2 tt 


so that 

h~ iw,4 + d^ V 1-df = 0. 

As expected in the limit £ — * 00 we obtain d\ ~ c -3i»/7£- 7 which is the limiting form of 
the solution obtained for ct ~ M ~ ’ , It is a simple matter to show from the above equation 
that dii is a monotonic decreasing function of l and in particular for £ — * 0 we find that 

Ji~rV" /3 . (3.20) 


The main implication of this limit is that the M~ * layer for the a ~ M * modes decreases 
in size like * when l\ -* 0. The next stage is when the M~* and M~ » layers merge 
which occurs when d\ = O(M) and a = 0(M~ 7 ). 

a = Q(M-») 

Here we expand c and a in the forms 

L 

c = co + . . . , a = — j + 

M 7 

Now there are just three regions, of depth M~ 7 ,M°, and M 7 , for us to consider. Inter- 
estingly the extent of the perturbation away from the wall has been reduced by a factor 
of M i - from the order M • scale of the case ct ~ M 7 . This reduction in depth is caused 
by (G - c) remaining 0(1) in this layer, and if we write 77 = M 7 Y, the dominant terms in 
the equation for the pressure in that layer give 

Pyy + T 2 (l — co) 2 p = 0 (3-21) 
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If co has a positive imaginary part, it follows that 


p = exp(-aY), 


with <7 = iL( 1 — co). Thus where Y —* 0 we have that ’ a. This means that 

the required solution for p in the 77 = 0(1) region is 


P= 1 + ~ TTrbi + »?] + (3.22) 

M 3 3 

Here the constant terms are small compared to M° but large compared to M~ » ; in prin- 
ciple they may be found by matching requirements. It remains to consider the small layer 
when 77 = M~ and p expands as 


# = 1 + -. + + 

P 2 is found to satisfy 

~ — ° P 2 - 9^ 2 (7 - 1)(1 + uo)(l - «o)(-(7 - 1)(1 + «o)(l - «o) - (wo - c) 2 ) 

tiQ CO Z Z 

so that 


p# = ( Uo _ c 0 ) 2 [ ( ^ 2 (7 ~ 1)(1 + wp)(l ~ wp) (7 - 1)(1 + w 0 )(l - up) 
*'« 2(uq-co) 2 5 


A match with the core solution is achieved if 


(wo - co) 2 )d£. 

(3.22) 



T 2 (7 ~ 1)(1 + wp)(l - up) 
4(wq - co) 2 


((7 - 1)(1 + wo)(l - wo) - 2 (uq - c 0 ) 2 K 


= — <7 = — iL(l — cq). 


(3.23) 

This is the required eigenrelation to determine the complex wavespeed co as a function 
of the scaled wavenumber L. A large L analysis of this equation shows that the limiting 
small £ solution of the a — 0(M~ 3) case is retrieved. 

In Figure (3.6) we show co as a function of L for 100 < L < 400 and the case of an 
insulated wall. We see that the maximum value of c^ for a two-dimensional mode occurs 


for L ~ 225 and a neutral mode exists for L ~ 110. We postpone until the final section a 
discussion of the implications of the results we have found above. 
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4 The Inviscid Instability Problem in Interactive Boundary Layers. 

We consider the hypersonic flow of a Sutherland law fluid past the semi infinite wall 
defined by the positive x axis. The oncoming flow has constant velocity at infinity and 
the leading edge of the wall is sharp. This results in the formation of a shock wave which 
acts as an upper bound for the flow disturbances in the sub-layers below, the uppermost 
of which is inviscid. The reader is referred to Figure (4.1) for a description of the different 
parts of the flow field. 

This region has been studied by for example, Bush (1966), and since viscosity is 
negligible here the use of a power- viscosity-temperature law does not alter the well-known 
governing equations for this region. Below the expansions and equations governing the 
flow are given. Apart from some minor differences, our formulation is essentially the same 
as that of, amongst others, Stewartson (1955, 1964) and Bush (1966) and so the reader is 
referred to those papers for more details. 

The flow quantities and coordinates have been non-dimensionalised using upstream 
flow quantities and L a typical streamwise lengthscale. We then scale the non-dimension- 
alised normal coordinate y appropriately 



(4.1) 


where M is the Mach Number and Y is taken to be 0(1). Note that as yet we do not 
scale the downstream variable x. As in previous studies we consider the two-dimensional 
problem and this enables the introduction of a streamfunction ip defined by 


pu = ip v , pv = -ip m 


(4.2a, b) 


Further we expand the velocity, temperature, density, pressure and streamfunction as 


U 1 

” = 1+ M2 + 


VI , 
V = -TJ + 

M 


T = Ti + 


P = Pi + 


Pi , 

p ~ M 2 *'"’ 


(4.3a, 6, c, d, e,/) 


and 


i> = 77^1 + 

M 


19 



If the above expansions are substituted into the inviscid Navier-Stokes equations written 
in Von-Mises coordinates we obtain, after some manipulation, the governing equations for 
this layer 

VU = , pi = E(ipi)p\. (4a, 6, c) 

The function E(ipi) is evaluated from the initial conditions which for these hyperbolic 
equations are given at the shock. Conventionally, we define the shock by Y = f(x) where 
/ is unknown at present; we note that 0i = f on the shock. 


The Strong Interaction Zone 

The solution of these equations (and the corresponding ones for the lower layers) can 
be investigated analytically in the large x— and small x— limits using expansion procedures, 
following Stewartson (1955, 1964), Brown and Stewartson (1975). The formulation below 
follows the latter paper closely. 

Let us consider the solution of (4.4a,b,c) for small x, i.e. close to the leading edge of 
the plate, where the shock is attached. Here we write 


/ = aix 3 ' 4 + • • ■ , pi — a?ft( Qx ~ 1/2 + • ■ • , 
v\ = aiu 1 (C)x _1/4 + ■ ■ ■ and Pi = Pi(C) H • (4.5 a,b,c,d) 


The similarity variable £ introduced in the above equations is defined by 


c = 


01 

a ix 3 / 4 


(4.6) 


The scalings (4.5b,c,d) are implied by the Rankine-Hugoniot relations, which relate 
flow quantities either side of the shock Y = f(x), and by (4.5a) which follows from matching 
normal velocities and pressure across the sub-layers beneath the shock; see Stewartson 
(1964) and references therein. 

The resulting equations for the leading order terms are 


with 


Pl«) = ^73^1(0 , VI + 3CUi = 4ft, 

_ 3C _ J 9(7 - ir 

Vl = 1 and 61 “ 8(7 + i)V+i’ 


(4.7a, b, c, d ) 


vi = 


3 

2?7+T) ’ 


Pi 


9 

8(7 + 1) 


Pi 


7+1 

7-1 


on C = 1 , 
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which are the conditions at the shock. We require their solution for £ < 1. In particular, 
for given 7 we can solve numerically for uj.(0) and Pi(0). 

Substitution of (4.5b,d) into (4c) and recalling the definition of C gi ves £(^1) ~ 
P\a\ 3 V>{ 2/3 as ipi -> 0. This implies that, as < -» 0, ~p x ~ C 2/37 and this produces a 
higher order correction in v\ if 7 > For a realistic gas 7 is larger than this value, 
with the result that there must be a viscous sublayer beneath the present layer. Bush 
(1966) found that this layer enables one to match the inviscid-layer solution on to the 
boundary layer solution if one used a power-viscosity/temperature law (/ j , ~ T w ,w < 1), 
rather than the linear law used by Stewartson and most other research studies in the field. 
Stewartson (1955) hoped that by using idealized assumed physical properties it would 
help to “understand the behaviour of more realistic fluids”. His assumption resulted in 
discontinuous derivatives, or “kinks” , in properties between the Inviscid Zone and the 
Boundary Layer, see Fig. 7.4, page 167 of Stewartson (1964). Despite the results of Bush 
(1966) (and a few subsequent authors) the linear viscosity law is still extensively used 
in hypersonic shock/boundary-layer research; it is argued that these ‘kinks’ are only a 
slight nuisance that can probably be explained away by deeper analysis and appeal to the 
argument of Stewartson concerning the assumption of ideal fluid properties. The analysis 
of Bush was extended to the linear case by Lee and Cheng (1969). They found that, by 
consideration of a second order boundary-layer correction, coupled to Bush s transitional 
layer analysis, the two layer structure of Stewartson gives the correct basis for the flow 
structure. 

Not surprisingly the extensive use of Chapman’s law in hypersonic calculations has 
encouraged researchers in stability theory to make the same rather severe approximation to 
the fluid viscosity; we believe that the calculations given in this paper are the first to take 
account of a realistic viscosity-temperature dependence in a first mode stability analysis. 
We note that our flow solutions in the layers beneath the shock are in full agreement with 
Bush, but he considered a general power law whereas we concentrate on n ~ T 1/2 , the 
leading order form for Sutherland’s law at high temperatures. For convenience we will 
use the notation of Stewartson and Brown (1975). We will next consider the viscous- 
adjustment-layer after discussing the third sub-layer, the viscous boundary layer. 

Viscous Boundary Layer 

This is adjacent to the surface of the plate and enables velocities to be reduced to 
zero on the surface. First, the full asymptotic formulation is described briefly, then the 
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asymptotic behaviour for small x is discussed in the next subsection. We assume that the 
boundary layer exists where ip ~ M -3 , ip being the stream-function defined earlier. We 
define the hypersonic parameter r by 


Re = rM 5 


(4.8) 


and take r to be order one. Note that the power of the Mach number is six in the definitions 
of the hypersonic parameters of those using a linear viscosity law. We write 


u = U i + 


Vi 

v = — + 

M 


P= Ml Pl ’ 


,, = w Rl + 


T = M 2 6\ and n = Mu i. 


(4.9 a,b,c,d,e,f) 


These scalings are standard apart from that for the viscosity n which follows immediately 
from the high temperature form of Sutherland’s law: 


H = (1 + C)T 1/2 . 


(4.10) 


We make the assumptions that the Prandtl number is unity and that the wall is 
insulating, this enables us to carry out a relatively simple analysis for the viscous solutions. 
The solutions capture the new structure of the inviscid modes due to the effect of using 
a power-viscosity law in the viscous sublayers. The relaxation of these assumptions does 
not substantially alter the conclusions obtained below. 

We again transfer to von Mises coordinates; using the streamfunction 'J' as an inde- 
pendent variable instead of Y. The energy and momentum equations are combined to 
yield the boundary-layer equations in these coordinates 


UiU u 


(7 ~ 1) (1 - U i 2 ) „ 7iWl 
2 7P1 lx r 


(Sr***) 


* 


and 

01 = ^(1 - V?), 

1/2 

together with /zi = B-y from (4.10), and with usual boundary conditions 


(4.11a, b) 


U\ = 0 on $ = 0, 

Ui -> 1 as $ -c 00. (4.12a, b) 


- » 
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Boundary-Layer Solution for the Strong Interaction Zone 

Now let x be small so that the hypersonic interaction parameter becomes large. We 
define the similarity variable <j) by 

* = Ml + OPlTOl 1 ' 2 ^) < 4 ' 13 > 

where a\ and Pi(0) were introduced when considering the inviscid region. If we now -write 

tfi = £(<£) + •••> ( 4 - 14 ) 


we obtain the similarity equation for the boundary-layer solution 


and 

g( 0) = 0 , g(oo) = 1 . (4.15a, b,c) 

This is a modified Faulkner- Skan equation, but the decay of g at large <f> is now 
algebraic rather than the usual exponential behaviour, 

a ~ 1 — —r + • • • , as <£->oo, ( 4 . 16 ) 

y (j t > 4 


and a > 0 can be determined analytically. 

Matching normal velocities across the sub-layers in the usual manner yields the leading 
order coefficient a\ in the small x expansion for the as yet unknown shock location 


« 1 


3(l + C) (' 2(7-m 1/2 /'7-l') 1/4 r 1-g 2 
4ui(0) ' r 7 p 1 ( 0 ) / ' 2 ' * o g 


(4.17) 


The Viscous Transition Layer 

Here the velocities and temperatures in the outer inviscid flow match with those in 
the inner boundary-layer flow. We give a brief account of the arguments used to deduce 
the position and properties of this region and again consider the strong interaction zone. 

Earlier it was noted that the inviscid density function p x ~ C 2/3t as C ^ can 

easily be deduced that u\ and T\ grow like inverse powers of C as C -♦ In fact 

T ~ Tl ~ x- 1/2 C 2/31 ~ ^ 1/2 ( ’ (4 ' 18) 
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as ip i — * 0 . 

Now let us consider the boundary-layer temperature as the top of the layer is ap- 
proached. We have written T = M 2 6i and found that 0i = ^(1 - U%). We know the 
large (p behaviour of U\ and it follows immediately that 


6\ ~ 


1 

4> 4 


tf 4 


X 

m^ipi ’ 


so that 


T = m 2 $ i ~ — £ 


7pim t 


as 


<f> — > oo. 


(4.19) 


Matching the limiting temperatures (4.18) (4.19) gives the position of the viscous transition 
layer, 


ipl ~ M~ x x( 3 ~ x ^* , 


X = 


97 
67 — 1 


(4.20a, b) 


This is in full agreement with Bush’s result. Note that the powers involved in the asymp- 
totics are functions of 7 , the ratio of specific heat capacities, in the new layer. Further 
investigation of the small ip\ limit of the inviscid solutions implies the following expansions 
for flow properties in the transition layer 


>lO 

u = l + M ••»-* «! + ••• , 
T = M •*- 1 T\ -(-••• , 


together with 


P = 

We define the scaled streamfunction 7 


IPX 
T ' 


(4.21a, b, c) 


ip\ = rjM A (= Mip). (4.22) 

These are then substituted into the Parabolised Navier-Stokes equations in Von- Mises 
co-ordinates to obtain the equations for u\ 

u l* _ TlPU 7Pl_g_f 1 

(1 + C) 7 pi r dr)\ j .7 dr)'' ' ° 

where, as we are taking the Prandtl number to be unity and the wall to be insulating, the 
corresponding energy equation can be integrated to yield 

T\ = (1 - 7 )« 1 - 
24 


(4.24) 


The pressure pi is independent of the normal co-ordinates, as usual in viscous regions 
governed by the boundary-layer equations, so the pressure is expanded 


pi = 


<*lPl(°) , 

r — + ••• 

x* 


in the two viscous layers, this being the limiting form as the inviscid layer is descended. 

We see that (from (4.23)) in this adjustment layer we still have the high temperature 
form of Sutherland’s law, in contrast with the shock-free adjustment layer of section 2 
where the full form is needed. Again, a similarity solution for the flow exists in this 
viscous sub-layer with similarity variable z defined by 


rj = /3x^z, 


(4.25 a) 


where 


p 2 _ (£ + l)7fliPl(0) 

T 


and 


Writing ui 
crucial region 


A = 


3 - A 
4 


(4.256, c) 


£ a- 2 G(z) i ea( i s to the governing equation for the flow solution in this 


(A - 2 )G - Q — ^ zG ' = i— + 
v ' 4 27 


G' 


V(l-7 )G 


(4.26) 


Note that from the definition of ui, (4.2 la), we expect that G < 0. We require that this 
solution for u\ matches to the corresponding solutions above and below, i.e. in the inviscid 
sub-layer and viscous boundary layer respectively, so we must investigate asymptotically 
the small and large z limiting forms of G. 

As z — ► 0 we find that G has the behaviour 

576 

G= -goz~ 4 + giz- v + go = ^ _ jp - (4.27a, 6) 


Thus the transition layer solution matches onto the large <(> form of g. The exponent u 
of z in the first correction term can be calculated analytically whilst the coefficient g\ 
must be determined from a numerical solution. In fact v satisfies a quadratic equation 
with coefficients functions of 7 . With 7 = 1.4 we find that u = 0.6267 ... so that we have 
relatively small correction terms. The value of go for this choice of gamma is 275.6 . . . and 
so G grows quickly as x — * 0. 
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We now investigate the behaviour of G for large z we find that 


G — > —Aqz~ + A\z 2 + • ■ • , 


(4.28) 


and so again we have a relatively small correction term. Note that G decays algebraically, in 
contrast to the rapid exponential decay of the Blasius and the “Modified Blasius” functions 
which arise in the shock-free, far-downstream cases (the former after employing the linear 
Chapman viscosity law and the latter from the use of the more realistic Sutherland’s law 
- see earlier). The second coefficient A\ is a function of gamma and can be determined in 
terms of .Ao; in fact 

A 8Aq (37 + 1) 

1 97(7 - 1)* (37 - 1) 

The constant Aq is not fixed by the asymptotics. Instead we must choose it so that the 
velocity solution matches with the appropriate limiting form of the inviscid solution. We 
obtain, after an elementary matching argument, and a little manipulation, the result 


Aq = (1 + C) 


1 _ j_ 2 — 


1 - ■ 


7 ]ej 

7-1 


(4.29) 


For the value 7 = 1.4 we find .Ao = 0.1751.... In our numerical calculations we took 
without loss of generality the hypersonic interaction parameter to be unity. Now that the 
base flow for this region has been found (at leading order) we can consider its stability 
characteristics; in particular we are interested here in the linear stability of inviscid modes 
concentrated (trapped) within this adjustment layer. These so-called vorticity-modes have 
been discussed earlier for other flows. 


The Vorticity Mode in the Strong Interaction Zone 

The scalings for these modes appear complicated but follow in a straight-forward 
manner after applying the usual vorticity mode arguments to the flow field discussed 
above. The modes of wavenumber k have wavelength comparable with the thickness of 
the transition region and a Rayleigh analysis suggests we require ~ T 2 k 2 . This can 
be seen, for example, from (2.12a) and (3.4) studied earlier. From (4.21b), and recalling 
the definition of A in (4.20b), we have T M 4a-6 , and also from (4.22), ij) ~ M x 1 . 
Note that with these scaling we are primarily concerned with obtaining the correct Mach 
number dependence; the small x dependence will be incorporated later. We deduce that 
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k ~ M >> 1) and thus this represents a short wavelength mode and so we introduce 
a fast x-scale 

X = M 7 - 3 x x. (4-30) 

The time-scale can also be similarly deduced from Rayleigh analysis, after noting that the 
vorticity modes propagate in a frame moving with uniform velocity and recalling the u\ 
sc alin g (4.21a) . We are lead to the introduction of a short time scale. 

r = M x ~ l t. (4.31) 


We note that these scales are compatible with our earlier analysis of the non-interactive 
case. At leading order we have the multiple-scales 

d z -+d t + M 7 ~ 3x d x , d t ^M x - 1 d T -M 7 - 3x d x . 


Note that non-parallel effects are 0(d x ) ~ 0(1) and are negligible in comparison with the 
direct growth effects of order k{u - 1) = 0{M X ~ X ) » 1; for our choice of 7 = 1-4 we 
find that A = 1.7027. The wavespeed c of the vorticity mode propagating in the frame 
moving with uniform (nondimensionalised) velocity u = 1 is related to its stationary-frame 
wavespeed c by 


c = 1 + 


M 4 (2-M 


(4.32) 


which follows immediately from (4.21a). Note that it represents a small perturbation about 
the velocity of the propagating frame. The remainder of the analysis foUows the classi- 
cal inviscid mode approach for formulating the pressure equation for the linear wave-like 
disturbance, now that the scales have been deduced. We assume that the flow previously 
discussed for this region is perturbed by an infinitesimal disturbance proportional to 


£ ik(X-cr) 


(4.33) 


containing all the X and r dependence of the disturbance. We note that the basic flow 
properties are independent of r and X allowing such an assumption. After a little manip- 
ulation we will obtain the so-called vorticity-made-Rayleigh equation for the amplitude of 


the disturbance pressure p v 


where we require p to decay to zero for large and small r\. We now return to the Strong 
Interaction Zone and investigate the asymptotic behaviour of (4.34) for small x. 

We recall the leading order small x forms 


PI 


= 2 ^ 21 , 


U l = X 


X 3 

x-2 


G. 


(4.35a, b , c) 


We need the small- a: dependence of c and k but need not worry about p as (4.34) is linear. 
We find that 

c = x x ~ 2 C, 

and 

r3 

k = ( 7 - < 4 ' 36 “' 4 > 
are the necessary behaviours. The former follows immediately from (4.35c) and the latter 
results from considering the coefficient of p in (4.34). We finally arrive at the renormalised 
problem for p; the vorticity-mode pressure- amplitude equation for the strong-interacticn 
zone 

2 G' o , 

Pxz - G _ £ p, - K*G*P - 0, (4.37) 

to be solved subject to p vanishing in the limits z — » 0 and z — » oo, so that, as in previous 
analyses, the disturbance is confined to the adjustment layer. We discuss our numerical 
solution of this eigenvalue problem and its limiting small wavenumber form in the next 
section. 
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5 The Solution of the Strong Interaction Vorticity Mode Equation. 

The leading order asymptotes of equation (4.37) are found to be 


and 


- 2 ( 

P~*P 0 z exp^-^ 3 -j as 

jl ( ZjAolC 
P -» PooZ" ex P y ~ 3 7 _ ~2 Z 


z — ► 0 

(37-2) \ 
37 ' 


asz — + oo 


(5.1a, b) 

where po and p & are constants and the numbers <70 and Aq arise from the corresponding 
asymptotes for G given earlier. Higher order terms in these expressions can be found 
analytically, and, in fact, are needed for accurate numerical solutions to this eigenvalue 
problem for C and /C, the scaled relative wavespeed and the scaled wavenumber of the 
inviscid vorticity mode, discussed after briefly paying attention to the neutral case. There 
C is real and takes the value of G evaluated at the generalized inflection point which of 


course occurs where 

GG lt = 2 Gl 


A numerical solution to (4.26) using a Runge-Kutta method shows that this occurs 
when f = 1.661432, where the new variable £ = tnz was introduced to ‘stretch’ the co- 
ordinate-scale for small z where G, and hence p, vary rapidly. The resulting neutral value 
of the wavespeed is then given by 


C = -0.633318. 


The corresponding real value of K is obtained from a numerical solution of (4.37) for all z 
with an appropriate treatment of the path of integration at the generalised inflection point 
(g.i.p.). The particular numerical procedure chosen to calculate this neutral wavenumber, 
and the complex eigenvalues C was to extend the line of integration into the complex- 
plane, taking a triangular contour around and below the g.i.p. The neutral value of K 
was calculated to be 


K = .477957 


Figure (5.1) shows the growth rate /CC* plotted against /C we find that the maximum 
growth rate ~ .060918 occurs at K ~ .156100 and that the growth rate goes to zero when 
the wavenumber goes to zero. In Figure (5.2) we show the eigenfunction of the most 
unstable mode, we note the exponential decay of the eigenfunction at both ends of the 
range of integration. 
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We shall now show how the small wavenumber structure of the strong interaction 
inviscid mode develops when the wavenumber tends to zero. In fact we simply give the 
structure corresponding to the first small wavenumber regime discussed in Section 3. Some 
discussion of the regimes encountered at even smaller wavenumbers will be given in the 
next section. 

The small-£ behaviour 

Numerical solutions of the pressure-amplitude equation, (4.37), indicate that C grows 
as K — » 0; exactly the same behaviour was found for the no shock problem studied in 
Section 3. However, in the latter, the leading order dependence on the wavenumber was 
indicated from the numerical solutions, whereas, in the present case, such guidance was 
not possible due to a complicated dependence of C on K. Instead, the equation was studied 
analytically for small K, and initial investigations suggested that 

1 4e J 

as k — » 0. Here the parameter e is a function of 7 and is defined by 

4e=l — (5.2) 

37 v ’ 

and we have replaced K, C by k, c respectively here and henceforth in this sub-section. 

4 

It follows from (5.2) that special attention is required for values of 7 close to 7 = — 

3 

and, as this coincides with the range of physically relevant values of the ratio of specific 

heat capacities, we solely consider the case of |e| << 1 in this section. We note that the 

4 

suggested form of c is defined and continuous at 7 = — by appealing to “L’Hopital’s Rule”. 

3 

Below we outline the asymptotic regions and corresponding solutions required to resolve 
the leading order behaviour of c as k — > 0. 

When z is very large, we have the WKB asymptote 

£, -L- ( 37^0 

P ~ z*y exp l kz »■» 

\ 37-2 

and we note that for smaller z this form breaks down in a ‘turning-point’ region corre- 
sponding to 

1 - 37 

z ~ k~ 3 -r-» >> 1. 

In this region we define 

y = k*^z, P = P 0 (y) + ■■■ (5.3a, b) 
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(5.4a, b ) 


and, using the large - z asymptotes 


for G, we find that Po satisfies 


Po;; - Aqv~ 4a Po = 0, 


a = 


37' 


This equation for Pq has an analytic solution involving the modified Bessel Function 


Jf , the required solution being 

P 0 = Do^K^ri) 


(5.5a) 


where 

77 = 4Ao^y 1 ', 1/ = 1 — 2a and ^ = — . (5.5 b>c,d) 

Here D 0 is an arbitrary constant. We can investigate the small y behaviour of Po by 
considering the limiting forms, for small arguments, of the modified Bessel functions. We 
find that 


Co = flo^2 ^[l + f^ 

Sin(/i7T) (1 — V) 


r( 2- M )(A^ 2 r , 

r(l-/x)(l-y) y 


(5.6) 


as y -> 0, the ordering of the second and third terms being dependent on the value of 7; 
we note that the corresponding powers of y (the indicia! roots of the associated Frobenius 
solution) become equal when 7 = — . Moreover, the coefficients of these terms are also 

O 

4 

singular when 7 has this value (note that 7 = - *-* n = 1). Here T is the Gamma function. 

That special attention is needed for values of 7 near four-thirds is now apparent and, 
as mentioned earlier, it is exactly this range of 7 that we are interested in, physically. 
Considering e = — — <x to be small, we see that 


Po; ~ 2A)Ao ( ~ 4g ~ ) + 1 + 0( c )] 


(5.7) 


as y -» 0. We note that P 0 ; (and, in fact, P 0 ) is defined, and moreover, continuous as 
e — 0. It is convenient to recast the last result in terms of z to ease matching to the next 
asymptotic region to be discussed. The result is that 


Pq z ~ 2D§Aftk 


k ( 


z 4e — 1 
4 F 


*• 

I? 1 - 3a 1 , 

)+ T - ) + l 

4 f 


(45) 


as z — * 0. 

The solution for P = Po(y) H , will continue to be the leading order behaviour of P 

until contributions from the ‘middle’-term (proportional to P,) become leading order. The 
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location of this next region of interest can be easily located, by analogy with the respective 
analysis of Section 3, as we expect c >> 1 as k -» 0. Here 2 ~ <r 1/4 , using the small-2 
asymptote of G, and we expand 

P = 2D 0 [l + k 2 6l /2 P(y) + •••], c = c6 2 + • • • (5.9a,. b) 

where 62 , = <^2(^1 f), is a small parameter. The P form is implied from the Pq solution for 
small y. The new, scaled normal variable, y is taken to be 0(1) and is defined by 

* = ^2 1/4 S (<< 1)- (5.10) 


Substituting into the pressure-amplitude equation, (4.37), yields the equation satisfied by 

P 

p— 1 8g ° p = 4 (x ,i\ 

vv V(go + cy 4 ) v y 8 ') 

where go arises from the small-2 asymptote of G. This can easily be integrated to give 


Py(V) = (c + 


^L) 2 ( [ v 

5 4 + go ) 2 



(5.12) 


with D\ an arbitrary constant to be determined by matching with the previous region. As 
V y 00 we see that Py — > Die 2 , thus P, — > .Dic 2 ^ 4 and matching with (5.8) gives 


k 2 Dic 2 6 7 J 4 = Atk^T' 


[( 


4e 

k 1 - 2a _ 1 
4 e 


M 


at leading order. We choose D\ = — Aqc 2 and this determines the size of c for small k 



(£7+77) 4e — 1 

— 4e ■ 


(5.13) 


4 

Taking the limiting value as e — ► 0, we have that when 7 = — 

= -lnfc 2 - 1 >> 1 (5.13) 

recalling that In x — > 00 as x > Oj the choice of Di <C 0 ensuring that 62 is real and we 

allow c to be complex. When e is small we see that 62 is large and has the form suggested 
at the start of this section. Thus 62, as given by (5.13), leads us to conclude that, for 
the physically realistic values of 7 around four-thirds, the scaled, relative wavespeed c 


32 


ii« in 11 1 1 


increases as the scaled wavenumber k decreases, this growth having logarithmic magnitude 
in contrast to the much faster algebraic growth found in the no-shock analysis presented 
in Section 3. We note that the corresponding frequency, proportional to ck, tends to zero 
as k decreases. The k- dependence of c, for k 0, can be deduced in a similar fashion 
for other values of and we note that, mathematically, the no-shock result of Section 3 

is the limiting case 7 — * 00 . 

AU that remains now is to determine c and this is easily deduced by introducing a 
further region where * ~ k 1 ' 3 « 1- The details are identical to those of Region IV in 
Section 3 ; again we find we need P,- -e 0 as S 0 , resulting in an expression for c 

r 9&v , 3. _ o 

J 0 (cf / 4 + so ) 2 C 2 

which can be solved for complex c; after an elementary contour integration and a little 
manipulation we obtain 


which corresponds to to an unstable mode of (4.37). 

Thus we have obtained the structure of the vorticity mode in the strong interaction 

region at smaU values of the wavenumber. We postpone any further discussion of our 
results until the following section. 


37t<7o 


(5.14) 
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6 Discussion 

We have investigated the instability of flat plate hypersonic boundary layers to the 
vorticity mode of instability. This inviscid mode is associated with the generalized inflec- 
tion point of the basic flow and is thought to be the most dangerous mode of instability 
of a high Mach number flow. W hen the mode is neutral the wave propagates downstream 
with the speed of the fluid at the generalized inflection point. At wavenumbers smaller 
than the neutral value the mode is unstable and the growth rate attains it’s maximum 
value at a finite value of the wavenumber. In the small wavenumber limit the growth rate 
approaches zero and for the non-interactive boundary layer at sufficiently small wavenum- 
ber the vorticity mode spreads out towards the lower boundary and reduces to an acoustic 
mode at a countable infinite set of wavenumbers. We believe that a similar process hap- 
pens in the strong interaction case since there the acoustic mode is correctly described by 
a quasi- parallel theory there. We did not pursue that calculation here because it would be 
essentially unchanged from that of Section 3 except that it would be made somewhat more 

complicated by the necessity of treating the case 7 = 4/3 as a special case in the strong 
interaction zone. 

We believe that the results we have presented in Section 4 are the first which show 
the effect of a leading edge shock on any form of hydrodynamic instability. Interestingly 
enough the shock does not have a direct influence on the vorticity mode; thus the main 
effect of the shock is to restructure the boundary layer in the leading edge region and 
thereby influence the susceptibility of the flow to inviscid disturbances. 

The vorticity mode eigenvalue problem was formulated in the interactive region along 
the plate at 0(1) values of x. However to the authors’ knowledge the basic flow in this 
regime has not yet been calculated; the numerical problem was set up by Bush (1966) but is 
sufficiently difficult to have remained unsolved. Thus we were unable to solve the eigenvalue 
problem in this regime and therefore choose to consider the strong interaction regime where 
a similarity solution for the basic state is available. An alternative to that limit would have 
been to consider the weak-interaction problem, Bush and Cross (1967), where a different 
similarity structure holds. We choose to concentrate on the strong interaction limit because 
the growth rates there are bigger and if the flow is indeed unstable there the stability of 
the flow further downstream is possibly of less relevance. 

Unfortunately we are unaware of any experimental observations or other theoretical 
work which we could compare with our results for the strong interaction regime. 
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In Section 3 we showed how the acoustic inviscid mode emerges from the small 
wavenumber description of the vorticity mode. Again it is not possible for us to com- 
pare our work with that of previous authors since it appears that the finite Mach number 
calculations available, mostly due to Mack, have either being carried out using a Chap- 
man viscosity law or a combination of Sutherland’s and Chapman’s law. In fact Mack s 
calculations were carried out using a combination of the different laws so as to efficiently 
model the viscosity-temperature structure of the fluid. The fact that the calculations of 
Cowley and Hall(1988) and Smith and Brown (1989) agree so well with Mack’s calculations 
suggests that over the part of the flow where instability took place Chapman’s law was 
being used; in the case of the vorticity mode this is clearly a bad approximation because 
the mode locates itself in the layer where the basic temperature field varies rapidly. 
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Figure (3.1) The dependence of the vorticity mode growth rate oic^ on the 
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Figure (3.3) The different regions which emerge in the small a limit. 
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Figure (3.6a) The real part of the wavespeed c 0 as a function of L for 
100 < L < 500 and the case of an insulated wall. 
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Figure (3.6b) The imaginary part of the wavespeed c 0 as a function of L 
for 100 < L < 500 and the case of an insulated wall 
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